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Delayed Enrichment by Unseen Galaxies: Explaining the 
Rapid Rise in IGM CIV Absorption from z = 6—5 
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ABSTRACT 

In the near future, measurements of metal absorption features in the intergalactic 
medium (IGM) will become an important constraint on models of the formation and 
evolution of the earliest galaxies, the properties of the first stars, and the reionization 
and enrichment of the IGM. The first measurement of a metal abundance in the IGM 
at a redshift approaching the epoch of reionization already offers intriguing hints. Be- 
tween z = 5.8 and 4.7 (a 0.3 Gyr interval only 1 Gyr after the big b ang) , the measured 
density of C iv absorbers in the IGM i ncreased by a factor of ~ 3.5 (|Rvan- Weber et all 
2009; iBecker. Rauch fc SargentJ[2009() . If these values prove to be accurate, they pose 
two puzzles: (1) The total amount of Civ at z = 5.8 implies too little star formation 
to reionize the IGM by z — 6 or to match the WMAP electron scattering optical 
depth (t). (2) The rapid growth from z « 6-5 is faster than the buildup of stellar 
mass or the increase in the star formation rate density over the same interval. We 
show that a delay of ~ 0.4-0.7 Gyr between the instantaneous production of ionizing 
photons and the later production of metal absorption features (added to the delay due 
to stellar lifetimes) can provide the full explanation for both puzzles. We calculate the 
delay in metal production due to finite stellar lifetimes alone and find that it is too 
short to explain the rapid C iv density increase. The additional delay could naturally 
be explained as the result of ~ 200 km/s outflows carrying carbon to distances of 
^100 kpc, the typical distance between galaxies and ClV absorbers in enrichment 
simulations, and the typical outflow or absorption region scale observed at z ss 2-3. 
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1 INTRODUCTION 

Most of the star formation taking place during the reioniza- 
tion epoch remains invisible to current observations, and a 
large fraction may remain invis ible even to the James Webb 
Space Telescope (JWST) (e.g. ISalvaterra. Ferrara fc Davall 
I20TQT ). This unobserved population of faint galaxies must 
have contributed most of the high-energy photons that 
reionized the intergalactic medium (IGM). The ionization 
state of the intergalactic medium is one important probe 
of this population. Other constraints are provided by the 
detection of luminous Lyman-break galaxi es (LBGs) and 
Lyman-q emitters (L AEs) at z > 6 (e.g. iBouwens et al.l 
l20ld : lTilvi et alfeoioh . Although the directly observed pop- 
ulation cannot account for all of the photons necessary 
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to reionize the IGM (and especially to match the WMAP 
Thompson-scattering option depth from f ree electrons r, 
iBolton fc Haehneltll2007l ; lOesch et al.ll2009l '), they constrain 
the high-luminosity end of the galaxy luminosity function 
(LF). 

The distribution (by element, in time, and in space) of 
the elements synthesised in these early generations of stars 
and subsequently expelled into the IGM (or incorporated 
into low-mass stars still observable in the local universe) will 
become an important source of information about the epoch 
of reionization in the near future. Currently the highest- 
redshift measurements of metal abundance in the IGM come 
from searches for C iv absorption features in quasar spect ra 
i|Rvan- Weber et al.ll2009l : IBecker. Rauch fc Sargendl2009l ). 

The uncertainties on the z > 5 C IV measurements are 
still large, but if the va lues prove to be accurate then they 
pose two puzzles. First, iRvan- Weber et alj (2009) find that 
the total amount of C iv at z = 5.8 implies too little star for- 
mation to reionize the IGM by z = 6. Second, IBecker et al.l 
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(2009) find that the fractional increase in C IV density is 
larger than either the buildup of stellar mass or the increase 
in the star formation rate density (SFRD) over the same 
interval. 

We show here that both puzzles can be solved by a delay 
between the production of ionizing photons and the enrich- 
ment of the IGM with the associated metals. In the next 
section (*J2J), we discuss current observations and models of 
C IV in the IGM in more detail. Section[3]presents our simple 
framework for modelling reionization and enrichment based 
on a galaxy luminosity function historyp In Section [4] we 
show that we can match reionization f i]4.1[l and enrichment 
constraints at a single redshift ( H4.2p with simple models and 
physically plausible parameters, then go on to show that a 
delay of 0.4-0.7 Gyr (in addition to the delay due to finite 
stellar lifetimes) can produce the observed rapid rise in C IV 
absorber density f i]4.3[) . In Section [5] we propose two expla- 
nations for the delay involving galactic outflows, and suggest 
observational tests of those explanations. In Section H3 we 
summarise our results and discuss how future observations 
and models will improve our understanding of early metal 
enrichment and reionization. 



2 EXISTING OBSERVATIONS AND MODELS 

Observations of absorption features in high-redshift quasar 
spectra are beginning to probe the IGM metallicity at red- 
shifts approaching z = 6 (arguably close to the final stages of 
reionization). Currently the highest-redshift measurements 
are of the Civ 1548.2 1550.8 A doublet redsh ifted into 
the near infrared fNIR) . iRvan- Weber et~afl (|2009h identified 
three (plus one tentative) C IV features in a search between 
z = 5.2 and 6.2 along lines of sight to 9 qua sars with a 
combi ned absorption distance of AX = 25.1rl iBecker et al.l 
(2009) performed a similar search in four sight lines (AX = 
11.3), finding no absorption features. This is particularly 
surprising since their obser vations were sensit i ve to even 
lower column densities than IRvan- Weber et al.l ([.2009), and 
therefore would have been expected to detect more C IV fea- 
tures per unit absorption distance if the column-density dis- 
tribution followed the declining power-law form found at 
lower redshift. 

Converting the absorption feature detections to the av- 
erage C IV density in the IGM (expressed as a fraction of 
the critical density), an d correcting for the i r com pleteness 
limits in column density, I Ryan- Weber et al.l (J2009J) found 



Apply ing similar corrections to observations bv lPettini et al.l 
(2003), they found 



fie iv {(z) = 5.76) = (5.0 ± 3.0) x 10" 



(1) 



1 The full compendium of code and ancillary files needed 
to reproduce the present paper is available from the 
first author. Cosmologial calculations were made with the 
CosmoloPy package. The EnrichPy package encapsulates 
our enrichment model. These resources are available at 
http : //www . astro . phys . ethz . ch/kr amer/ 

http : //roban . github . com/CosmoloPy/ and 
http://roban.github.com/EnrichPy/ 

2 AX is defined so that objects with constant comoving density 
and physical cross section have cons tant density per unit AX 
[|Tvsonlil988l ; IRvan- Weber et al.ll2009l) . 



ficiv ((*) = 4 -69) = (17 ± 6) x 10" 



(2) 



The errors are still large on these measurements. If the 
z = 5.76 value is revised toward the upper end of the al- 
lowed range, then the puzzle of the rapid evolution in C IV 
will be greatly reduced (see Figure U in Q . However we 
will proceed here under the assumption that the central val- 
ues are essentially correct. We therefore seek to resolve the 
puzzle of the rapid C IV growth. 

The C IV density in the IGM is appr oximately constant 
at fi c iv ~ 2 x 10~ 8 from z = 2-4 (|Rvan- Weber et al.1 
120091 ; ISongailal 120011 ). a surprising result given that this 
is a period of intense star formation. Models of IGM en- 
richment have successfully explained this as the result 
of a decreasing fraction of carbon in the triply-ionized 
state, offsetting the concurre n t rise in total carbon density 
dOppenheimer fc Davd l2006l : iDave fc Qppenheimerl 120071 ; 
lOppenheimer fc Davell200Sl ). 

Models of enrichment and carbon ionization at 2 > 
4 have been able to produce a rise in ficiv be- 
tween z = 6 and 5 consistent with the observa- 
tions (|Oppenheimer. Dave fc Finlatorl 120091 ; ICen fc Chisaril 
20101) . at least w i thin their large uncertainties. 



Both 



Qppenheimer et all (|2009l ) and ICen fc Chisaril IJ2010D find 



that the total amount of carbon in the IGM increases by a 
factor of 2.5 to 3 in this interval, while the fraction of car- 
bon in the triply-ionized state only increases by a factor of 
^ 1-25. 

These simulations do not self-consistently model reion- 
ization, so they are unable to directly elucidate the relation- 
ship between reionization and enrichment, and therefore un- 
able to provide a satisfying solution to the first puzzle posed 
above, which is our primary concern here. We are therefore 
interested in investigating the connections between reioniza- 
tion, enrichment, and the population of galaxies responsible 
for both. 



3 METHODS 

3.1 Connecting Reionization and the Galaxy 
Luminosity Function 

In order to solve the two puzzles posed by the rapid fic iv 
rise between z — 6 and 5 (JEJ, we need to model the lumi- 
nosity function of galaxies, and from that model calculate 
the total ionizing emissivity as a function of redshift. We as- 
sume the LF is a Schechter function. Observations suggest 
that <f) t (the density normalisation) and a (the faint-end 
slope of the LF) are a pproximately constant from z — 4 to 6 
IjBouwens et al.1120071 ). Neither are constrain ed as precisely 
at higher redshift, but iBouwens et al.l (2010) find that the 
observed LF at z — 7 and 8 is consistent with constant a 
and 0*, though the maximum-likelihood values move to- 
ward steeper faint-end slopes (a = —1.94 ± 0.24 at z ~ 7 
and —2.00 ± 0.33 at Z ~ 8). We therefore parameterise the 
evolution of the LF by varying onl y M t with reds h ift. A t 
Z < 9.0 we use the M, values from IBouwens et al.l {2008), 
linearly interpolated as a function of redshift. At z > 9 we 
assume M* is linear in redshift with slope /3m*. The param- 
eters a and /3m* therefore control the extrapolation of the 
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LF to lower luminosities and higher redshifts than have yet 
been probed by observations. 

To convert the observed LF to an ionizing photon 
emi ssivity (or rate den s ity, p hotons -1 s _1 Mpc 3 ), we use 
the iBolton fc Haehnerq (|2007m spectral energy distribution 
(SED) and an escape fraction / esC 7- With this SED, the ratio 
of ionizing photon production rate to UV luminosity is 



N 



L(1500 A) 



: 8.4 x 10 



24 photons s 
erg s _1 Hz" 



(3) 



though we assume only a fraction / eS c7 of these escape into 
the IGM0 

With these factors we convert the integrated LF his- 
tory into an ionizing emissivity history, and then into the 
ionized fraction of the IGM x(z). We ignore twice-ionized 
helium, assume that once-ionized helium has the same num- 
ber fraction as hydrogen, and include recombinations of hy- 
drogen. Recombinations are calculated with clumping factor 
C = (n^n) / (nuii) 2 — 4 (see H4.ll for discussion of varying 
C), gas temperature 10 4 K (giving case B recom bination rate 
a B = 2.6 x 10~ 13 cm 3 s' MHui fe Gnedirj[l997h and assum- 
ing all ionized gas is contained in fully ionized bubbles. We 
integrate the LF down to Mab{UV) = — 13.04, equivalent 
to a star formati on rate of 0.01 mQ/yr (IKennicuttl 1 19981 ; 
lOeschet al.ll2009l ). We use the "W MAP7 + BAO + HO" 
mean cosmological parameters from iKomatsu et al.l (|20101 ) 
throughout this paper. The values of / OSC 7, /3m*, and a are 
discussed in §3] 



3.2 Connecting Reionization and Carbon 
Production 

We define f x civ as the ratio between the total rate (with 
no delay) of C IV production and the total rate of ionizing 
photon production, so that: 



fie 



t(z) = X to t*\{z)f x CIV, 



(4) 



where fie ivinst(z) would be the rate of C IV density increase 
with no delay in emission, and itotai(z) is the ratio of the 
total rate density of ionizing photon production to the total 
number density of hydrogen and helium atoms. 

The fxciv ratio depends on a number of factors: 



fxCIV = 



/civ/ei 



r 7 z 



(Xc/Z)Q 



bar yon ■ 



(•">) 



The fraction of all metal m ass in carbon is Xc/Z = 
0.178 (the solar value from lAsplund. Grevesse fc Sauvaj 
|2005|).fibaryon = 0.0456 is the fraction of the cr itical density 
contributed by baryons (jKomatsu et al.l 120101 ). The terms 
in square brackets (the fraction of carbon in Civ, /civ; the 
fraction of metals that escape galaxies, f csc z', and the ra- 
tio of metal nucleon to ionizing photon production, r 7 z) are 
highly uncertain. We use the representati ve values of /civ = 
0.5 (the maximum t heoretical value, see iRvan- Weber et al.l 
120091 ; ISongailal 120011 ) and f CS cZ = 0.2 (equal to our fiducial 



3 Note that 7V/L(1500 A) is sensitive to the initial mass function 
and mctallicity of the stellar population, though we ignore this de- 
pendence here. 7V/L(1500 A) is completely degenerate with / eS c7 
in our formalism. 



/csc7), and emphasise that it is only the product of these 
uncertain factors that matters. 

The ratio of ionizing photons to metal nucleons pro- 
duced by a stellar population can be expressed as 

E n 



r-,z 



V 



E a 



(6) 



where 77 is the ratio of total ionizing photon energy to total 
rest-mass energy of the metals produced in a stellar popu- 
lation, E p is the rest-mass energ y of a proton, and -E a v g is 
mean energy of ionizing photons (Schaerer 2002). For a stel- 
lar population with Z = 1 /5O/?0 . ISchaererl (2002) calculate 
77 = 0.014 and £ avg = 21.95 eV, yielding f xC iv = 1.4x 10~ 9 . 
Note that this 77 value does not include yields from low- 
and intermediate-mass stars (LIMS, m < 8ttiq), stellar 
wind mass loss, or Type I SN contributions, all of which 
would decrease 77 (and increase the enrichment to ioniza- 
tion ratio fxciv)- We will also explore scenarios using the 
solar-metallicity (Z = Zq) values of 77 = 0.0036 and 
-Eavg = 20.84 eV (calculated with mass loss and SN Ibc, 
but still without LIMS), resulting in f x civ = 5.0 x 10~ 9 . 

Changing the M* vs. z slope (/3m*) or the faint -end 
slope of the LF (a) changes both the ionization and enrich- 
ment histories. Changing / OS c7 only affects the ionization 
history. Changing f x civ only affects the enrichment history. 



3.3 The Carbon Delay Distribution Due to Stellar 
Lifetimes 

Carbon is not produced instantly upon formation of a pop- 
ulation of stars, unlike, for our purposes, ionizing photons. 
Instead carbon is ejected primarily after the main-sequence 
lifetime of a star is over, which for low and intermediate mass 
stars (m < 8Mq, tm e > 3x 10 7 years) becomes a signifi- 
cant fraction of the relevant timescales (e.g. the 3 x 10 s year 

interval from z = 6 to 5). 

Neglecting this delay (as in lCen fe C hisari 201CJ) is often 
justified with the statement that Type II supernova (SNII) 
yields from short-lived, high-mass stars dominate carbon 
production at high redshift. However, even if most most 
carbon is synthesised in SNelJj , a substantial fraction of 
that carbon is incorporated into low- and intermediate-mass 
stars (LIMS, m < 8mQ) b efore being blown into the IG M, 
according to the models of lOppenheimer fc Pavel (J2008J). In 
other words, an important fraction of carbon produced in 
a galaxy gets locked up in LIMS and is only returned to 
the gas phase (and made available for ejection in a galactic 
outflow) during the asymptotic giant branch (AGB) phase, 
after the main-sequence lifetime has elapsed. Therefore, the 
lifetimes of lower-mass stars impose a delay on the ejection 
of some of the carbon into the IGM, and AGB star ejection 
of carbon cannot be neglected in calculating the timing of 
IGM enrichment. 



4 The fraction of carbon contributed by low and intermediate 
mass stars depends sensitively on the amount of "hot bottom 
burning" (HBB) that takes place on the asymptotic giant branch 
(AGB), since HBB can destroy carbon and even result in a net loss 
of 12 C over the lifetime of a star. The amount of HBB as a func- 
tion of stellar mass is still highly uncertain, iVentura &: Marigd 
I2010T) so it is un clear whether low— o r high— mass stars dominate 
12 C production JRomano et alj|2010l . and references therein). 



© 2010 RAS, MNRAS 000, 000-000 



4 Kramer, Haiman & Madau 




0.00 
10 



10 7 10° 

delay t (years) 



10 3 




0, 8.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 



delay t (years) 



le9 



Figure 1. Delay functions for carbon production. Differential 
(top panel) and cumulative (bottom panel), delays are calculated 
with Chabrier (thin curve) and Kroupa (thick curve) IMFs. The 
differential delay is shown on a logarithmic time scale, while the 
cumulative delay is shown on a linear time scale. Note the rapid 
emission from supernovae (~ 10 7 years) and the longer tail from 
AGB stars (~ a few xlO 8 years). 



Since we are trying to model the connection between 
ionizing emission associated with star formation and the 
eventual enrichment of the IGM with carbon produced by 
the same stars, it is important that we take this delay into 
account. Indeed, such a delay is inevitable, and the original 
motivation of this paper was to assess whether this delay 
might help explain the steep observed evolution of the C IV 
abundance in the IGM, as discussed in Sjl] 

The cumulative delay function fc(t) is the fraction of 
carbon emitted by stars with lifetimes tm e < t. Assuming 
all carbon is ejected at the end of a star's main sequence 
lifetime, the fraction of carbon ejected in the time interval 
t ± (dt/2) after star formation is 



dfc , ■, _ dfc dm 
dt ^ ' ' ~ dm dt' 



(7) 



where dm/dt is the inverse of the derivative of the lifetime 
function (tiif e (m)). The fraction of carbon produced by stars 
of mass m ± (dm/2) is 



dm 



(m) oc Mc{m)ij>(m), 



(8) 



where Mc(m) is the carbon mass ejected by a star of mass 
m, and <j>(m) is the initial mass function (IMF) of stars by 



number. We normalise this function to give J max dfc/dt — 
1. Here t min = 3.24 x 10 6 years is the lifetime of a 1OOM 
star, and i max = 1-52 x 10 9 years is the cosmic time between 
the starting and ending points of our simulations, z = 100 
and z = 4.1. With our chosen IMFs and yields (see below), 
an additional < 9% of carbon would emerge at t > t max . 

Convolving the fie ivinst (i) curve with 3L£ gives us the 
delayed C IV curve 

fie IVdclay(i) = fie IVinst * ( —j— J , (9) 

whe re * represents convolution over the time coordinate. 

iRomano et al.l (|2005l . l2010h have quantified in detail the 
impact of uncertainties in the IMF, stellar lifetimes, and 
stellar yields on Galactic chemical evolution models. The 
greatest uncertainties are in the yield calc ulations, which 
vary c onsiderably from author to author. As Romano et al. 
(20 1G) point out, there is no consistent set of yields covering 
the whole range of mass and metallicity and including all 
of the physical effects relevant to either galactic or cosmic 
chemical evolution models, and essentially no suitable cal- 
culations have been performed for m m 6 - 8Mq . For conve- 
nience, we use a set of carbon ( 12 C) yield values provided by 



IGavilan. Buell fc Mollal (|2005l l 5 l , contai ning their own orig- 
inal c alculations for m = 1-8 Mq , and IWooslev fc Weaver! 
(1995) values for m = 8 - 1OOM . We use the yields for 
metallicity Z = 0.02 because the variation in calculated car- 
bon yield with metallicity (at least above some threshold) is 
smaller than the overall uncertainty in yields. 

The final complication with the use of yield tables is 
the distinction between the total mass of an element ejected 
by a star and the net yie ld of new atoms synthesised in the 
star (IGavilan et al.l l2005) . In a self-consistent chemical evo- 
lution model that tracks the metallicity of the star-forming 
environment, the yield of new elements is the relevant quan- 
tity. However, we are only tracking carbon abundance in the 
IGM, while stars are forming directly out of gas in the in- 
terstellar medium (ISM) , so we use the total ejected mass of 
carbon to calculate the delay. Note that this is independent 
of the calculation of the total amount of carbon produced 
(see i )3.2[) . This is a good approximation if the ISM reaches a 
stable metallicity quickly, and the composition of the galac- 
tic outflow is representative of the total mass currently being 
ejected from stars (both via SN and AGB mass loss). 

Theoretical calculations of stellar lifetimes generally 
agree fairly well for m > IMq. The dependence on metal- 
licity is quite weak. The larger uncertainties at m < 
IMq are irrelevant here since the corresponding lifetime 
of tnie > 10 Gyr is longer than the age of the universe 
at z > 0.3. We adopt the lif etime function of Kodama as 
given in lRomano et al.l (j.2005). 

The stellar initial mass function is another important 
source of uncertainty. We can characterise the impact of the 
IMF on the delay function by dividing the carbon emis- 
sion into a prompt component and a delayed component. 
We are concerned here with evolution on a timescale of 
~ 10 8 years. Therefore carbon emission that occurs faster 
than 10 7 years after star f ormation is rel a tively prompt. Of 
the IMFs discussed in the lRomano et al.l ((2005J) review, the 



Available on VizieR: http://vizier.u-strasbg.fr 
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iKroupa, Tout fc Gilmorel ||l993h and IChabrierl l|2003h IMFs 
produce the most extreme values for the cumulative de- 
lay function at 10 years, /c(10 yr). With a Kroupa IMF 
/c(10 7 yr) = 0.23, while fc(W 7 yr) = 0.44 for a Chabrier 
IMF. We therefore calculate all of our models with both of 
these IMF^j in order to demonstrate quantitatively the im- 
pact of the IMF uncertainty on our results, and to suggest 
qualitatively the impact that different sets of yield values 
might have. Figure [TJ shows the differential and cumulative 
delay functions with Kroupa and Chabrier IMFs. These fig- 
ures illustrate that roughly half of the carbon ejection occurs 
essentially instantly, whereas the remaining half is spread 
over the lifetime (~ lGyr) of intermediate-mass stars. 

Extremely low metallicities may result in dramatically 
different IMFs and yields from those assumed here. We do 
not consider this metal-free, or Population-Ill, mode of star 
formation in calculating the delay function, as metal-free 
stars are thought to make up only a small fraction of the 
stars formed before z = 6, even if their formation continues 
at a low rates to late times (see, e.g. iRollinde et al.l 120091 ; 
ISalvaterra et alJl2010h . 



4 RESULTS 

The primary constraints on the epoch of reionization avail- 
able today are the WMAP measurement of r (the optical 
dept h to Thompson scatt ering from free electrons in the 
IGM. lKomatsu et alj|2010f ) and the evo lving L yman-a opac- 
ity of the IGM at z ~ 6 (though see iMesin gcr 2010 for a 
discussion of the complicated, model-dependent interpreta- 
tion of this evolution). In i]4.1l we discuss the agreement 
between our luminosity function histories and these reion- 
ization constrains. Then in sections 14.21 and 14.31 using the 
formalism outlined above, we return to the puzzles posed in 
the introduction. 



4.1 Matching Reionization Constraints 

In this section, we compare our luminosity function his- 
tories to available constraints on the reionization of the 
IGM. For our fiducial luminosity function history, we adopt 
a Schechter luminosity function with fixed <j>* = 1.1 x 
10~ 3 Mpc~ 3 a n d a, and use the observed M* values from 
iBouwens et al.l (|2008h (linearly interpolated) from z — 3.8 
to 9.0. We extrapolated M, linearly in z above z — 9 with 
slope Pm*- The extrapolation t o lower luminositie s is con- 
trolled by the faint-end slope, a. lOesch et al.l (120090 suggest 
that M* is approximately lin ear in z at high redshift, with 
a slope of /3m* = 0.36 ± 0.18. IBouwens etail (120081 ) suggest 



a = —1.74 at high redshift. We use f e , 



0.2. 



Figure [5] shows the ionization and enrichment histories 
as functions of redshift. The lower (blue) set of curves corre- 
sponds to the fiducial parameters described above. The top 
panel shows the ionized fraction in the IGM, along with the 
optical depth integrated from redshift to z. The arrows 



6 Note that we are only varying the IMF in the calculation of the 
delay function. In principal N / L(1500 A) and r 1 z also depend 
on the IMF and metallicity of the stellar population, but we treat 
each of these calculations independently. 
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Figure 2. Ionization and enrichment histories of the IGM. In each 
panel, the lower (blue) set of curves corresponds to the fiducial 
LF history, and the upper (red) set to the enhanced LF chosen to 
match WMAPt = 0.087 and z re ion ~ 6.6. Top panel: ionized frac- 
tion x(z) (solid) and optical depth r(< z) (dotted). The arrows 
on the right show the total value of the optical depth (integrated 
to z = 100). Middle panel: Civ density in the IGM with (solid) 
and without (dashed) a delay function. The thin solid curve uses 
the Chabrier— IMF delay, while the thick curve uses the Kroupa- 
IMF delay. Observed values are indicated as black data points 
with error bars. Horizontal error bars indicate the interval over 
which the density is averaged, and vertical error bars indicate the 
uncertainty. Bottom panel: Enrichment in the same models as in 
the middle panel, except with f x civ increased by a factor of 3.7. 
This increase is equivalent to the difference between 1/50-solar 
and solar— metallicity values for the ratio of ionizing photons to 
metal nuclcons. 
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on the right show the tot al value of t h e opt ical depth (in- 
tegrated to z = 100). As Oesch et al.l 1)20091 ) have pointed 
out (see also lBolton fc Haehne lt 2007), combining this lumi- 
nosity function evolution with reasonable ionizing photon es- 
cape fraction (/ OS c7 = 0.2) and IGM clumping factor (C = 4) 
values, yields an insufficient emissivity to either complete 
reionization by z ~ 6 or match the WMAP constraint on 
the electron scatter ing optical depth (r = 0.087 ± 0.014 
iKomatsu et al.ll2010l ). With the fiducial LF history, the IGM 
is only 20% ionized by z = 6, and the optical depth is 
t = 0.036. 

In order to match the WMAP r value and complete 
reionization at z > 6, we increase / CSC7 iV/I/(1500 A) by a 
factor of 3, flatten the M* slope to /3a/* = 0.09, and set 
the faint end slope to the steeper value of a — —1.95. The 
higher / CSC7 iV/Z/(1500 A) could be explained by a higher 
escape fraction, or a higher N/L(1500 A) due to a lower 
metallicity or mor e top-heavy IMF of the stellar popula- 
tion. For instance, ICharvl (J2008I ) finds a factor of ~ 2-3 in- 
crease in JV/L(1500 A) when the metallicity falls from OAZq 
to 0.02Zq. Note that such changes would affect the delay 
function and r 7 z values as well. The changes to /3m* and a 
make t he enhan ced L F history resemble the recent results by 
iBouwens et al.l (|2010l ). who found a brighter M* and steeper 
q at z = 7 and 8 than were suggested by earlier results. We 
integrate the LF down to Mab(UV) = — 13.04, equivalent 
to a star formati on rate of 0.01 m,Q/yr (|Kennicuttl 1 19981 ; 
lOesch et al.ll2009l ). Adjusting any of these parameters alone 
cannot match both constraints, and adjusting them simulta- 
neously allows us to use more plausible values. More impor- 
tant than the exact parameter values is the resulting ionizing 
emissivity history. This is obviously not a unique solution, 
but we present this enhanced LF history as a plausible ex- 
ample of one that matches current observational constraints 
much better than the fiducial extrapolation of the observed 
LBG luminosity function. The upper (red) sets of curves 
in Figure [2] correspond to this enhanced LF history. The 
top panel shows that the optical depth has been increased 
to near the WMAP value, and the IGM is fully ionized by 
2 = 6.6. Table [T] summarises the parameters of each LF 
history. 

We use a constant clumping factor C — 4 in these calcu- 
lations, though our basic conclusion that the enhanced LF 
history is consistent with existing reionization constraints 
is not particularly sensitive to changes in this assumption. 
For instance, using a higher constant C — 6 results in 
t — 0.067, Zrcion = 6.21 for the enhanced LF history. The 
clumping factor should actually be lower at higher redshift, 
however. ICharvl 1)20081 ) has derived the clumping factor as 
a function of redshift for the relevant gas (ionized gas out- 
side of ion izing-photon source halos) from simulations by 
iTrac fc Cenl 1)20071 ). Using their clumping factor history (es- 
timated from their Figure 2a), results in a larger r = 0.082 
(because the clumping factor is lower, C < 4, at early times, 
Z > 10.5) and a later jz rc i on = 5.7 (because C > 8 at 
2 < 7 ). This optical depth is quite close to the WMAP7 
value. IPawlik. Schave fc van Scherpenzeell ()2009al lbh found 
that reheating of the IGM results in an even lower clump- 
ing factor history, which would require less enhancement in 
the LF history (something between our "fiducial" and "en- 
hanced" LFs) to match reionization constraints. 



Table 1. Parameters of our luminosity function (LF) histories. 



a &M* fesc-f 



fiducial -1.74 0.36 0.2 


1.2 


0.036 


enhanced -1.95 0.09 0.6 


6.6 


0.075 


WMAP7 (Komatsu et al. 


2010): 


0.087 ±0.014 



Parameters: 

a is the faint-end slope of the Schechter luminosity function. 

/3m* is the slope of M* as a function of z. 

fcecj is the escape fraction of ionizing photons from galaxies. Note 

that it is completely degenerate with the SED slope, which we fix 

at iV/L(1500 A) = 8.4 X 10 24 photons s _1 /(erg s" 1 Hz" 1 ). 

Results: 

^rcion i s the redshift at which the ionized fraction x(z) = 1. 

T is the optical depth due to free electrons. 

4.2 Matching CIV Abundance 

The bottom two panels of Figure [2] display the C IV mass 
density as a fraction of the critical density Clcw(z). The 
dashed lines assume instantaneous ejection, the thin solid 
lines are convolved with the Chabrier-IMF-based delay 
function, and the thick solid lines with the longer Kroupa- 
IMF delay. 

In order to show the effect of the large uncertainty in 
the value of f x civ, in the middle panel we use the 1/50 Zq 
value for f x civ, while in the bottom panel we use the 
solar-metallicity value (see i]3.2[) . The solar-metallicity stel- 
lar population produces fewer ionizing photons per metal 
nucleon synthesised, so for a fixed LF history it produces 
higher fie iv values. Both of the puzzles discussed in Sjl]are 
evident in the mid dle panel, which effectivel y uses the same 
fxciv assumed by I Ryan- Weber et al.l (| 20091 ) . 

First we ca n see the conflict suggested by 
iRvan- Weber et al.l 1)20091 ) between reionization con- 
straints and the low z = 5.8 Qc iv value in the middle panel 
of Figure [2] The enrichment history calculated from the LF 
history that matches reionization constraints (upper/red 
dashed curve) produces too much Civ at z = 5.8lJ This 
overproduction at a single redshift, in and of itself, is not too 
troubling. Either a slight decrease in the highly-uncertain 
fxciv or a slight decrease in the star formation rate density 
(SFRD) could lower flc iv sufficiently to agree with the 
2 = 5.8 measurement. Also, the stellar-lifetime delay, 
especially with the less top-heavy Kroupa IMF, brings 
f2c iv down to close to the observed value. 

The second puzzle — the rapid buildu p of C iv from 
2 = 5.8 to 4.7 noted bv lBecker et al.l (120091 ) — is also evi- 
dent in Figure [2] In the middle panel, again, the enhanced 
LF curve with no delay is close to matching the z = 4.7 
observation, but slightly over-predicts the earlier z = 5.8 
point. Similarly, the fiducial LF curve with no delay only 
slightly under-predicts the z — 5.8 value, but is much lower 
than the later z = 4.7 observation. The predicted evolution 



7 Again, we are ignoring the large errors on the observations for 
the sake of exploring their consequences should they prove to be 
accurate. If we take the observational errors into account, then 
we can see that even the upper/red dashed curve is marginally 
consistent with the observations. 
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of fie iv using either LF history is too slow to match the 
observations. 

We were motivated to calculate the stellar-lifetime- 
based enrichment delay by the idea that such a delay might 
help to explain to rapid rise in Civ. In fact, as Figure [2] 
clearly illustrates, we find that stellar lifetimes contribute 
little to the solution of this puzzle. In the next subsection, 
we explore the effect of longer delays. 

4.3 The Rapid Rise in CIV 

Since f x civ is so uncertain, it is useful to find a quantity 
independent of f x civ to compare with the observations. For 
a given combination of an ionizing emissivity history and a 
delay function, the fractional increase in fie iv over a spe- 
cific redshift interval is fixed and does not depend on f x civ 
(as long as f x civ is independent of redshift). The observed 
fractional increase is 



2.4. 



(10) 



Aficiv _ ficiv (4.7) 
ficiv ~ ficiv (5.8) 

To compare this number with our theoretical curves, 
we average fie iv over the same intervals used to determine 
the observed values. With no delay, the fractional increase 
is Afi c iv/fic iv = 0.79 for the fiducial LF history, and 0.70 
for the enhanced LF history. These values are far smaller 
than the observed increase. 

We expect a delay in enrichment to make the fractional 
increase in Civ larger, because (if ficiv evolves slower than 
exponentially) the fractional rate of increase at earlier times 
must be higher. To see this, consider that the fractional in- 
crease is roughly 



Afi 



civ 



fit 



-(t)At, 



<> ■ «) •■'-■ C 11 ) 

i.lC XV "CIV 

where At ~ 0.3 Gyr is the time interval from z — 5.8-4.7 
and fie iv is the derivative of fie iv with respect to cosmic 
time. If fie iv is a power-law, proportional to t n (this is a 
reasonable approximation at the cosmic times we are con- 
sidering, and it is also a conservative one, in the sense that 
structure formation is increasingly more rapid at higher red- 
shifts), then 



fi 



civ 



fi 



CIV 



(*) 



t 



(12) 



Therefore, the fractional growth rate decreases with time. 
If we had a delta-function delay, so that fie ivdoiay (t) = 
fie ivinst(i— ideiay), then the delay would boost the fractional 
increase by a factor of 



t 



(13) 



fie IVdclay/fic IVdclay 
fic IVinst/fic IVinst « ~~ *delay 

We need a factor of t/(t — tdelay) = 3-3.5 to boost the frac- 
tional C iv increase from Afic iv/fic iv = 0.7-0.8 to 2.4. We 
can conclude that a delta-function delay of ~ 0.7 Gyr should 
boost the growth sufficiently to match the slope inferred 
from the two observations. Our stellar-lifetime-based de- 
lays, however, are too short to provide the necessary boost. 
The mean delay with the Chabrier IMF is 0.16 Gyr. With 
the Kroupa IMF, the mean delay is only slightly longer, 0.25 
Gyr. Furthermore, and even more problematically, the de- 
lay is not a delta function. A large fraction (roughly half; see 



Table 2. Fractional increase in Civ density for each luminosity 
function history and delay function combination. 



LF History 


Delay Function 


/xC/v/10- 9 


Aficrv/^civ 




(1) 


(2) 


(3) 


observed Clc 


IV evolution 


— 


2.4 


enhanced 


no delay 


1.5 


0.7 


enhanced 


Chabrier 


1.9 


0.7 


enhanced 


Kroupa 


2.3 


0.8 


enhanced 


0.5 Gyr delta 


4.5 


1.1 


enhanced 


Cha + 0.5 Gyr 


6.2 


1.3 


enhanced 


Krp + 0.5 Gyr 


7.8 


1.5 


enhanced 


0.6 Gyr delta 


5.9 


1.1 


enhanced 


Cha + 0.6 Gyr 


8.4 


1.8 


enhanced 


Krp + 0.6 Gyr 


10.7 


2.1 


enhanced 


0.7 Gyr delta 


7.9 


2.1 


enhanced 


Cha + 0.7 Gyr 


11.8 


2.6 


enhanced 


Krp + 0.7 Gyr 


15.7 


3.1 


fiducial 


no delay 


3.0 


0.8 


fiducial 


Chabrier 


3.9 


0.9 


fiducial 


Kroupa 


4.7 


1.0 


fiducial 


0.5 Gyr delta 


10.8 


2.3 


fiducial 


Cha + 0.5 Gyr 


16.3 


2.7 


fiducial 


Krp + 0.5 Gyr 


21.7 


3.2 


fiducial 


0.6 Gyr delta 


15.6 


4.0 


fiducial 


Cha + 0.6 Gyr 


24.7 


4.8 


fiducial 


Krp + 0.6 Gyr 


34.6 


5.4 


fiducial 


0.7 Gyr delta 


25.4 


8.2 


fiducial 


Cha + 0.7 Gyr 


42.7 


9.5 


fiducial 


Krp + 0.7 Gyr 


62.5 


10.5 



(1) Figure[3]shows the corresponding enrichment history for each 
row in this table (except the 0.6 Gyr delays). 

(2) The f x civ value given is the one needed to match Cc IV a t 
z = 4.7. Compare to our estimates of f x civ = 1-4 X 10 — 9 for 
Z = 1/50 Z , or 5.0 x 10" 9 for Z = Z@. 

(3) AS7civ/^CIV is independent of f x civ- 



Figure [TJ of the carbon is ejected promptly, which dilutes 
the boosting effect. 

Figure [3] shows enrichment histories with f x civ ad- 
justed to fit ficiv(-z = 4.7) to the observed value. This 
allows us to assess the effect of the delay on the slope by 
comparing the earlier evolution of the ficiv curve to the 
observed value at z — 5.8. If C IV absorber production is 
too slow, Afic iv/fic iv will be too small, and we will over- 
predict the earlier measurement. Table [5] gives the f x civ 
and Afic iv/fic iv values for each combination of LF his- 
tory and delay function. In Figure [4] we plot Afic iv /ficiv 
as a function of the delay. 

The stellar-lifetime delays alone are insufficient to boost 
the Civ growth to the observed rate, whereas assuming a 
0.7 Gyr delay instead brings the enhanced LF history into 
agreement with the observations (top panel of Figure [3}. 
Since the stellar-lifetime delay is inevitable, even if a longer 
delay is also in effect, we next add the 0.7 Gyr delay to 
the stellar-lifetime-based delay (middle panel of Figure [3} . 
When the stellar-lifetime delays are added, the fractional 
growth is increased by 25-50% (for the Chabrier and Kroupa 
delays, respectively). In the bottom panel of Figure[3]we add 
a shorter 0.5 Gyr delay. Table[2]also includes an intermediate 
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Figure 3. Enrichment histories with the stellar carbon— to- 
ionizing— photon production ratio, f x ciVi adjusted to match the 
z = 4.7 observation. The dashed curves assume instantaneous 
production of C IV. Top panel: comparing stellar-lifetime-based 
delays and a 0.7 Gyr delta function delay. Middle panel: a 0.7 Gyr 
delay has been added to the stellar-lifetime delays. Bottom panel: 
a 0.5 Gyr delay has been added to stellar-lifetime delays. 
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Figure 4. Fractional increase in C IV density as a function of the 
additional delay t. Note the log scale of the vertical axis. The 
bottom three curves show the fractional increase in C IV density 
from z = 5.8 to 4.7 for the enhanced luminosity function his- 
tory with no stellar-lifetime delay, the Chabrier-IMF-based de- 
lay, and the Kroupa-IMF-bascd delay (bottom to top). The top 
three curves show the corresponding increases with the fiducial 
LF history. The delay t is in addition to any stellar-lifetime de- 
lay. The solid horizontal line indicates the observed value of the 
increase (2.4), while the dotted lines show the (lcr) confidence 
limits on Qciv(z = 5.8). The curves with the additional stellar- 
lifetime delay (top two in each set) match the observed value at 
t = 0.4 to 0.7 Gyr, depending on the LF and IMF. 



0.6 Gyr additional delay. Figure fallows an easy comparison 
between the predicted and observed values of Af2c rv/flc iv- 

The stellar-lifetime delays have an effect on both the 
length of the additional delay and the f x civ value needed 
to match the observations. For instance, with the enhanced 
LF history and a delta function delay of 0.7 Gyr, the frac- 
tional increase is 2.1 and f x civ — 7.9 x 10 -9 . With the 
Kroupa-IMF delay included, a shorter additional delay of 
0.6 Gyr produces the same fractional increase, but then 
fxciv must be increased by 35%. Therefore, future mod- 
els designed to study the relationships between reionization, 
enrichment, and galaxy formation must include the finite 
stellar lifetimes in order to draw precise quantitative con- 
clusions. 

While some of the f x civ values in Table[2]are physically 
plausible, some may be unphysically high, indicating that 
certain delay function and LF history combinations are not 
reasonable candidates for explaining the observed increase 
in C IV density. For instance, the fiducial LF with a delta- 
function delay of 0.7 Gyr requires the rather high value of 
fxCiv = 2.5 x 10~ 8 . On the other hand, f xC iv = 7.9 x 10~ 9 
(for the enhanced LF history with 0.7 Gyr delay function 
delay) could plausibly be explained as the result of a solar- 
metallicity r-,z (ionizing photon to metal nucleon ratio) and 
fescz = 0.3 (instead of 0.2). 

So far in this section, we have ignored the observational 
uncertainties on fie iv- We are less concerned about the error 
at Z — 4.7, since Qc ivis approximately constant at z < 4.7. 
Therefore in Figure f4] we indicate the range of the frac- 
tional increase corresponding to the lcr confidence interval 
on f2civ at Z — 5.8 (dotted lines). The current large ob- 
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servational errors allow a wide range of values for the delay 
(e.g. la limits of ~ 0.1-0.7 Gyr for the fiducial LF), but 
Figure |4] makes it clear that the current constraints exclude 
zero delay at more than the la level, and that a combina- 
tion of tighter constraints on fie rv{z = 5.8) and the galaxy 
luminosity function at z > 6 has the potential to place in- 
teresting constraints on the delay. 

In Figure [4] we have shown that, without any change 
in the C IV ionization correction, a ~ 0.4-0.7 Gyr delay be- 
tween the production of ionizing photons and C IV absorp- 
tion features can explain the rapid incr ease in C IV density 
betw e en z = 5.8 and 4.7. Simulations l|Oppenheinier et al.l 
12009 ; ICen fe Chisaril 120101 ) suggest that the triply-ionized 
fraction of carbon may increase by a factor of up to 1.25 in 
this interval, which would help to explain the observed rise, 
but would leave a factor of > 2.7 increase (Afic/fic = 1-7) 
in the total carbon content of the IGM. Figure|4]tells us that 
this growth still requires a delay of ~ 0.6 Gyr with the en- 
hanced LF history, or ~ 0.4 Gyr for the fiducial LF. In the 
next section we will explore possible physical mechanisms 
for such delays. 



5 DISCUSSION 

Since stellar lifetimes were too short to provide the required 
delay, what mechanism could explain a longer ~ 0.4-0.7 Gyr 
timescale? We show in this section that changing the stellar 
initial mass function is not a viable explanation, but that 
galactic outflow timescales correspond nicely to the required 
delay. We then make testable predictions for two different 
outflow-driven delay scenarios. 



5.1 Failed Explanations for the Delay: the IMF 

One way to produce a longer mean delay in carbon produc- 
tion would be to change the stellar IMF to increase the pro- 
portion of long-lived, low-mass stars|j To produce the long 
delay times found above requires an IMF that is radically 
bottom-heavy at high redshift. Even with a steep power- 
law IMF of <j>(m) oc m" 3 ' 5 , we find Aficiv/ficiv < 1-1. 
Therefore we also explored using a truncated Salpeter IMF 
((f>(m) ex m -2,35 ) with different maximum masses, and found 
that a maximum stellar mass of 2.44 rriQ (with the fidu- 
cial LF) or 2.09 rriQ (with the enhanced LF) is required 
to produce a Aficiv/ficiv of 2.4. These stellar masses 
correspond to lifetimes of 0.4 and 0.6 Gyr, respectively, 
in agreement with the delay times that we found above. 
However, these IMFs would produce no ionizing radiation 
to accomplish reionization, and would not match the con- 
straints on the SEDs of high redshift galaxies. Further- 
more, a bottom-heavy IMF is in opposition to the trend 
expected for stars forming from m etal-poor gas (but see 
lOmukai. Schneider fe Haiman 2008, who propose that when 
the metallicity exceeds a critical value of ~ 10~ 5 Zq, dense 



8 Though th is is outside the stellar mas s range we are con- 
cerned with, Ivan Dokkum fe Conrovi 2010 have found evidence 
for <f>(m) oc m -3 at m < 1 jtoq in massive local early-type galax- 
ies. 



clusters of low-mass stars may form in the nuclei of second- 
generation galaxies). Therefore changes to the IMF seem 
totally unable to explain the rapid fie iv evolution. 



5.2 Explanations for the Delay: Galactic Outflows 

In order to explain a long delay between star formation and 
C IV absorber production, we may posit that carbon ejected 
from a galaxy must reach a characteristic distance before 
producing a C IV absorption feature visible in current data 
sets. In this case, the relevant timescale for a distance d and 
velocity v would be 



t = 0.5 Gyr 



d 200 km s" 



(14) 



100 kpc v 

The characteristic distance from the source galaxy at which 
C IV is observed can be affected both by the ionization state 
of the carbon as a function of distance, and the filling factor 
of observable C IV absorbers. If the volume filling factor of 
the absorbers is too low, absorption systems will become too 
rare for detection in the limited set of sight lines currently 
available. 

While the exact relationship between galaxies and ab- 
sorbers is still highly uncertain (and contested), observa- 
tions at z = 2-3 s e em c onsistent with this outflow sce- 
nario. ISteidel et al.l (|2010l ) studied absorption associated 
with Lyman-break galaxies, both along the line of sight to 
the galaxies and along the line of sight to background galax- 
ies at impact parameters from 3-125 kpc. They found that 
the absorber equivalent width in their composite spectra 
declines slowly with impact parameter (Wo(CTV) oc b~ a2 ) 
up to b ~ 80 kpc, after which it quickly drops below their 
detection limit (Wo ~ .1 A). The absorbers detected by 
iRvan- Weber et all J2009D at z = 5.8 range from W = 0.06 
to ~ 0.7 A, which would place them at impact parameters 
from b > 60 to b > 100 kpc. They fit the Wo versus 6 pro- 
file with a model of a spherically-symmetric outflow with a 
covering fraction by absorbing clouds of f C ov{r) oc r -0 ' 23 . 
The slow decline in covering fraction with distance indicates 
that the absorbing clouds must be expanding as they move 
away from the source galaxy (otherwise the covering frac- 
tion would go like r -2 ). For a constant expansion velocity 
and conserved cloud number Re oc r 1 '' 2 when / cov oc r -7 . 
The inferred cloud radius therefore increases as Re oc r ' 9 . 
The fact that the C IV equivalent width declines more slowly 
than other ions (for which / cov oc r -0 ' 6 ) may indicate that 
ionization effects are also serving to increase the C iv/C ra- 
tio with distance from the source galaxy. Therefore, these 
observations suggest that both an increasing filling factor of 
absorbers and an increasing triply-ionized fraction of car- 
bon with distance from the source galaxy could delay the 
app earance of C iv abso rption in the IGM. 

ISteidel et al.l (|20100 also use their direct line of sight 
(impact parameter 6 = 0) absorption profiles to constrain 
models of the outflow velocity. They reproduce the observed 
profiles with an accelerating velocity profile (higher veloci- 
ties farther from the source galaxy) reaching maximum ve- 
locities of ~ 800 km s _1 . Given the ~ 100 kpc distance 
corresponding to a b sorbe rs with the equivalent width of the 
IRvan- Weber et all (J2009J) sample, this gives a time scale of 
> 0.1 Gyr. With the higher sensitivity of quasar spectra (but 
using the absorber-absorber correlation function rather than 
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direct detection of the associated galaxies), iMartin et al.l 
(2010) find a larger size for absorption regions of ~ 150 kpc, 
though they also suggest that the size may be constant in 
comoving coordinates, which brings the physical size back 
dow n to 85-100 kpc at z = 6-5. 

IMartin et ail i|2010l ) and lTvtler et all <J2009ft . using the 
distortion between the line-of-sight (redshift space) and 
transverse correlation functions of quasar spectrum ab- 
sorbers constrain the peculiar velocity of the average ab- 
sorber to < 200 km s~ . However, these analy ses compared 
the m ean or central velocities of the absorbers. ISteidel et al.l 
(J201CJ) showed that the absorption profiles are quite pre- 
cisely centred on the redshift of the source galaxy (at least 
at low impact parameters and averaged over the angular 
size of the background galaxy), with the width of the profile 
extending to high velocities. Therefore these results would 
constrain only the peculiar velocities of the host galaxies, 
not the o utflow vel o cities. 

Both ISongailal (|2006T l and iFox et all (|2007l ) measured 
the velocity widths of C IV absorbers, and found them to 
be correlated with absorber column density. In the former 
sample, absorbers with column densities comparable to the 
z = 5.8 sample have widths (at one-tenth maximum), of 
20-300 km s" 1 , while in the later (a study of DLA- and 
sub-DLA-associated absorption), all of the absorbers with 
IogJV(CTV) > 13.6 have widths of > 100 km s _1 , and most 
are 200 km s _1 -400 km s" 1 . 

Numerical hydrodynamic simulations also seem to sug- 
gest t hat both ionization and fi lling-factor effects are at 
work. lOppenheimer et al.l (|2009T ) find that C IV absorbers 
at z — 6 tend to have left their source galaxy 0.1-0.5 Gyr 
earlier, and lie within 10-50 kpc of a galaxy in their sim- 
ulation (giving average speeds of ~ 30-300 km/s). They 
suggest that the ionizing radiation from the source galaxy 
large ly determines the dist ance at which carbon is seen as 
C IV. ICen fc Chisaril (|2010T ) found absorbers at similar dis- 
tances (~ 70kpc). In their simulations, collisional ionization 
is important to the C IV ionization balance, and they sug- 
gest that at least some C IV absorbers are produced in shocks 
for med in galactic outflows. 

lOppenheimer et al.l (|2009h also find that, as metals 
travel farther from the source galaxy, they enrich less dense 
regions of the IGM and produce weaker C IV absorption sys- 
tems. Therefore (neglecting the ionization effects for the mo- 
ment) absorbers near a galaxy will have a higher column 
density, but a smaller filling factor (and will therefore be de- 
tected more rarely), while absorbers farther from the galaxy 
will have a lower column densities, but larger filling factors. 

Whatever is determining the distance scale, a ~ 0.4- 
0.7 Gyr delay could naturally be explained as the result of 
~ 200 km/s outflows carrying carbon to distances of ~ 80- 
140 kpc before it is seen in Civ absorption. 

5.3 Predictions 

We have shown that a delay between the production of ion- 
izing photons and of (currently-) detectable C IV absorbers 
in the IGM can solve the puzzles posed by the rapid rise 
in Oc iv between z ~ 6 and 5, and have suggested that the 
delay could be due to either (or to a combination) of two 
mechanisms associated with galactic outflows. Each mecha- 
nism leads to distinct observable predictions. 



If the filling-factor evolution dominates the delay, it 
leads to an interesting prediction: if the C IV has indeed al- 
ready been ejected from galaxies, and the reason it is not 
yet seen at the highest redshifts is that is has not yet spread 
to occupy a detectable filling factor, then future observa- 
tions, with significantly larger effective lines-of-sight, should 
uncover much of this hitherto hidden carbon in rare, high- 
column density Civ absorbers (at low impact parameters 
to unseen galaxies) |f| This model would predict that, once 
the rare, high-column-density end of the column density dis- 
tribution has been probed, the total integrated ficrv curve 
will be a scaled version of the dashed curves in Figure[31 and 
the f^c ii curve will be a scaled version of the f^c iv curve. 
Together, these observations would confirm that the amount 
of C IV present in the IGM tracks the cosmic star-formation 
history, but the rise in the filling factor of C iv systems is 
driven on a longer timescale, determined by the finite (rela- 
tively low) speed of the carbon-transporting winds or out- 
flows. The evolution of the C iv column-density distribution 
and correlation function will provide constraints on outflow 
and enrichment models. 

On the other hand if ionization effects (either local 
or universal) determine the delay, then ongoing and future 
measurements of the density of C II should reveal that much 
of the missing carbon is hidden in that ionization state. Sim- 
ilarly, observations of Sill, Si ill, and Siiv absorbers should 
reveal a changing ionization balance in that element. 

These are clear and feasible ways to discriminate obser- 
vationally between ionization-driven evolution and outflow- 
filling-factor-driven evolution. In either case, a careful ac- 
counting of the total amount of carbon in the IGM (sum- 
ming over ionization states and the absorber column density 
distribution) will reveal whether it tracks the cosmic star- 
formation history. Any departure from a constant propor- 
tionality either indicates a problem with measurements of 
the star-formation history, or a changing efficiency of carbon 
production, both of which would be of considerable interest. 



6 CONCLUSIONS 

We have shown that a ~ 0.4-0.7 Gyr delay between the pro- 
duction of ionizing photons and C iv absorption features can 
explain the rapid evolution of the C iv density in the IGM 
between z = 5.8 and 4.7. No change in the ratio of carbon 
to ionizing photon production (the ionizing efficiency) or in 
the universal ionization correction for C iv is required. This 
delay has a natural physical explanation, namely the need 
to transport carbon a certain distance into the IGM before 
it is seen in C iv absorption. The distance scale would be de- 
termined by a combination of the need to enrich a sufficient 
volume of the IGM to be detectable in a limited number of 
quasar sight lines, and possible ionization effects that opti- 
mise the C iv fraction at a certain distance from the source 
galaxy. An outflow of 200 km s" 1 would carry material to a 
distance of 80-140 kpc on these timescales. 

The shorter delay due to finite stellar lifetimes cannot 



9 The lack of detection of weak C IV absorbers in the lBecker et al.l 
(2009) observations is evidence against the alternative hypothesis 
that the carbon is hidden in systems too weak to detect with most 
current spectra. 
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provide the full explanation for the rapid evolution (even 
with a steep stellar initial mass function), but must be in- 
cluded in future models in order to properly understand the 
relationships between galaxy formation, reionization, and 
enrichment of the IGM. 

Future measurements of metal abundance in the 
IGM will provide important constraints on these rela- 
tionships. Perhaps the most important improvement will 
come from probing more quasar lines of sight, both to 
reduce the uncertainty in the density of absorbers in 
the currently-detected c olumn-density range (10 13,5 -10 14,5 ; 
lltvan- Weber et alJ 12009). and to search for the rare, high- 
column-density absorbers that we expect to exist close 
to galaxies. This search does not req uire extraordinarily 
high resolution or signal-to-noise ratio l|Oppenheimer et al.l 
|2009j), just the identification of more high-z quasars and the 
taking of their near-IR absorption spectra (though probing 
the weak end of the column-density distribution will require 
higher quality spectra). The Mul ti Unit Spectroscopic Ex- 
plorer (MUSEi lfBacon et alj|200q ). an integral field unit for 
the VLT, will be ideal for directly exploring the relation- 
ship between galaxies and IGM absorbers because it will 
be able to simultaneously obtain spectroscopic redshifts for 
large numbers of galaxies (up to z ~ 6.6) near a quasar 
line of sight, which can then be compared with the redshift 
distribution of absorbers seen in the quasar spectrum. 

Observations of multiple ionization stages of several ele- 
ments, using the sensitivity and wavelength coverage of new 
and future near-IR spectrographs like X-shooter, will break 
the degeneracy of the factors entering into f x civ, since some 
factors should be (approximately) independent of element 
(such as f C scZ and r y z), while the fraction of carbon in the 
C IV stage (/civ) can be constrained from measurements 
of Cn absorption and of the io nization balance in other 
species (|Rvan- Weber et al.l |2009f ) . Detailed physical mod- 
els will be crucial in interpreting these future observations . 
Such work is alread y underway IjOppenheimer et al.l 12003 ; 
ICen fc Chisaril 120101 ) , but needs to be improved to model 
both reionization and enrichment self-consistently. 
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